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Abstract 

We report on the astrometric observations of the 22 GHz H2O masers in the high mass star- forming re- 
gion G5.89-0.39 with VERA (VLBI Exploration of Radio Astrometry). Newly derived distance of 1.28t°;o8 
kpc is the most precise and significantly nearer than previous values. We revised physical parameters and 
reconsidered nature of G5. 89-0. 39 based on the new distance as follows. (1) The ionizing star of the ultra 
compact (UC) Hll region is a late 0-typc (08 - 8.5) zero age main sequence (ZAMS) star, consistent with 
previously established limits based on its infrared spectral line emission. (2) Crescent-like maser alignment 
at the position of the O type ZAMS star may trace accretion disk (or its remnant), which suggests that 
the star is still young and before complete evaporation of circumstellar materials. (3) Although the revised 
mass for the east- west outflow has been reduced, it still quite large (100 Mq) which indicates that a signifi- 
cant fraction of the mass is entrained material and that the dynamical age significantly underestimates the 
actual outflow age. Our newly-derived distance emphasizes that G5. 89-0. 39 is one of the nearest targets 
to investigate ongoing high-mass star formation and evolution in a compact cluster containing a young 
0-type star. 

Key words: Astrometry — ISM:masers (H2O) — ISM individual (G5.89-0.39, W28A2) — ISM:HlI 
regions — ISM:jets and outflows 



1. Introduction 

A systematic and consistent scenario of high mass star 
formation has not been constructed yet. In spite of enor- 
mous and intense works, there are many hypotheses and 
unsolved issues (e.g., Zinnecker & Yorke 2007 and refer- 
ence therein). But, several recent theoretical studies have 
suggested that high mass star formation can be achieved 
via mass accretion (e.g., Krumholz et al. 2009). This hy- 
pothesis seems to be consistent with observational signa- 
tures of massive accretion disk and torus around high mass 
protostellar object (HMPO) or protocluster (e.g., Beuther 
et al. 2009). Upcoming Atacama Large Millimeter / sub- 
millimeter Array (ALMA) will be able to resolve such a 
circumstellar structure enough and provide us more quan- 
titative information about a specific accretion mechanism 
onto an individual HMPO. 

Precise distance of a source is essential to quantitative 
discussions. However, most of the high mass star-forming 
regions, which will be targets for ALMA, are even located 
on the inner Galactic plane, where source distances of- 
ten contain significant uncertainty. Direct distance mea- 
surements by astrometric observations are quite impor- 
tant for such regions, in especial, highly accurate VLBI 
(Very Long Baseline Interferometry) astrometry is the 



only technique suitable for deeply embedded high mass 
star-forming regions where optical measurements are al- 
most impossible. Hachisuka et al. (2006) has been first 
demonstrated great performance of VLBI astrometry us- 
ing H2O maser in W3(0H) region, and since then, several 
studies have achieved 10 micro-arcsecond (/zas) accuracy 
for northern star-forming region (e.g., Moellcnbrock et al. 
2009; Sato et al. 2010b). In this paper, we report on an 
annual parallax measurement of H2 O masers in the high 
mass star-forming region G5.89-0.39 with VERA (VLBI 
Exploration of Radio Astrometry; Kobayashi et al. 2008). 

G5.89-0.39 (also known as W28A2) is one of the most 
famous, shell type ultra compact (UC) Hll region (e.g.. 
Wood & ChurchwcU 1989). The O-type ionizing star has 
been detected as a near-infrared (NIR) point source by 
Feldt et al. (2003) inside the shell (hereafter Feldt's star). 
Acord, Churchwell & Wood (1998) (hereafter ACW98) 
have directly measured dynamic angular expansion of the 
radio shell. Observed supersonic expansion and short dy- 
namical age (600 yr) indicate that this small UCHlI region 
is just after the birth. G5. 89-0. 39 is also known to be a 
host of an extremely massive outflow which is centered 
on the shell (e.g., Acord, Walmsley & Churchwell 1997). 
The whole part of the shell is completely inside the outflow 
extent (e.g., Watson et al. 2007). This also gives further 
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support on a remarkable youth of the UCHlI region. 

Previously reported distances for G5. 89-0. 39 vary over a 
wide range (1.9 - 3.8 kpc; Hunter et al. 2008, and reference 
therein). Almost all of them arc measured through kine- 
matic distance method, but at the Galactic longitude of 
5°. 89, this method intrinsically contains large systematic 
error of kpc order. Although ACW98 tried to estimate the 
distance from the shell expansion without any Galactic ro- 
tation models, they still have adopted several assumptions 
for their modeling of the data. In this point of view, our 
direct distance measurement is very important to confirm 
the physical parameters of G5. 89-0. 39. 

2. VERA Observations 

VERA observations of the H2O masers {Jk^k^ = 616- 
523) at 22.23508 GHz associated with G5.89''-0'39 have 
been carried out at 9 epochs between 2007 November and 
2009 May. A summary of all observations is listed in table 
1, which contains observing dates, synthesized beam sizes 
in milli-arcsecond (mas), beam position angles (PA) east 
of north, spectral resolutions, typical noise levels and dy- 
namic ranges of synthesized images and brief comments, 
if any, about system noise temperatures (Tsys)- 

Each observation was made in VERA's dual-beam mode 
in which targeted maser source and phase calibrator 
(or position-reference source) were observed simultane- 
ously (Kawaguchi et al. 2000; Honma et al. 2003). The 
real-time, instrumental phase difference between the two 
beams was measured for a calibration using the arti- 
ficial noise sources during the observations (Honma et 
al. 2008a). We chose J1755-2232 (a2ooo=17''55™26?2848, 
i52ooo=-22°32'10"61656; Petrov et al. 2005) as a paired 
calibrator. This source is separated from G5. 89-0. 39 by 
1°.92 at a position angle of -140° east of north. The flux 
density of J1755-2232 is about 150 mJy. No significant 
structure is seen in this point-like calibrator. A bright 
calibrator, NRAO530 (=J1733-1304; Ma et al. 1998), was 
also scanned every 120 minutes as a delay and bandpass 
calibrator. Each observation was made for about 6 hours, 
but total scan time for the targeted source pair was only 
about 2.5 hours because we observed not only G5. 89-0. 39- 
J 1755-2232 pair but also another maser-calibrator pair 
alternately. We will describe another source pair in forth 
coming paper. 

Left-handed circular polarized signals were quantized at 
2-bit sampling and filtered with the VERA digital filter 
unit (Iguchi et al. 2005), after that, data were recorded 
onto magnetic tapes at a data rate of 1024 Mbps. There 
were 16 IF channels with 16 MHz band width where one 
IF was assigned to the maser lines and other 15 IF of to- 
tal 240 MHz were assigned to J1755-2232 and NRAO530. 
The data correlation was performed with the Mitaka FX 
correlator (Chikada et al. 1991). Correlated data were 
divided into 512 and 64 spectral channels for the maser 
and calibrators, respectively. For the maser lines, the full 
16 MHz data were used only in the fourth epoch, and in 
other cases, wc used 8 MHz which is covering whole maser 
emissions to achieve sufficient velocity resolution (0.21 km 



s-i). 

The system noise temperatures were depended on 
weather conditions in each station and elevations. Typical 
Tsys value at a averaged elevation angle was varied 200 
to 500 K for each station in a case of normal weather. 
Sometimes it exceeded 1000 K in a bad case for Ogasawara 
and Ishigaki station. There was very high Tgys value (~ 
5000 K) at Iriki station in the first epoch, and hence, 
we performed phase-referenced imaging using other three 
stations in this epoch. Phase-referenced images of the tar- 
geted maser source were successfully obtained in 8 epochs 
without the final epoch, where significant phase fluctua- 
tions were still remained after the fringe fitting for J 1755- 
2232, and then, all maser spots were completely defocused. 

3. Data reduction 

Data reduction was carried out using the NRAO 
Astronomical Imaging Processing System (AIPS) package. 
Amplitude and bandpass calibrations were made for the 
targeted maser and J1755-2232 independently. We cali- 
brated clock parameters for each station using the residual 
delay of NRAO530. The tropospheric zenith delay offset 
was also calibrated by the modified delay-tracking data 
which were calculated based on the actual measurements 
of the atmospheric zenith delay with the global position- 
ing system (GPS) at each station (Honma et al. 2008b). 

There were two different paths of the analysis depend- 
ing on a purpose. The one was phase-referencing (or 2- 
beam) analysis for a measurement of the annual parallax 
and another was single-beam analysis which was suitable 
for searching all maser spots and studying internal mo- 
tions of them. Here, the term 'maser spot' indicates a 
maser component seen in a single velocity channel, and 
by contrast, the term 'maser feature' represents physical 
gas clump which is consist of several maser spots which 
are detected in successive velocity channels and closely lo- 
cated each other. We define it in the same way as Motogi 
et al. (2008). 

In the single beam analysis, we simply performed fringe 
fitting and self-calibration using the brightest maser spot 
in the feature 01 (sec section 4.3). We then searched 
for all maser spots by a wide-field mapping under a 7-a 
detection limit. Total explored area was 5x5 arcsec^ of 
each channel map centered on the ionizing star in G5.89- 
0.39. Annual aberration was also corrected for several 
masers which were significantly distant (^ 3") from the 
phase-referenced maser spot in this analysis,. 

In the phase-referencing analysis, fringe fitting and self- 
calibration for J 1755-2232 were done, and then, obtained 
delay, rate and phase solutions were applied to the tar- 
get visibilities following the correction of measured in- 
strumental phase difference between the two beams (see 
above) . The modified delay-tracking mentioned above was 
re-calculated with respect to each maser feature for ac- 
curate position measurement. New delay-tracking center 
was always located within 10 mas from relevant maser 
feature. 

The coherence in phase-referenced images were still sig- 
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Date 


Synthesized Beam 
(mas X mas) 


PA 

(°) 


Resolution 
(km s-i) 


1-cr Noise" 
(Jy beam^^) 


DN„'' / DNs 
(ct) 


Comments 

for Tgys 


2007/11/5 


1.92 X 0.71 


-24.33 


0.21 


0.82 


9/32 


^ 5000 K at IR= 


2008/1/12 


1.88 X 1.02 


-16.23 


0.21 


0.35 


11 / 166 


- 1000 K at IS 


2008/3/14 


2.15 X 0.88 


-20.15 


0.21 


0.22 


15 / 402 




2008/5/7 


2.32 X 0.88 


-24.48 


0.42 


0.63 


14 / 530 


500 - 2000 K at OG 


2008/7/1 


2.14 X 0.81 


-23.04 


0.21 


0.42 


8 / 128 




2008/11/11 


2.62 X 0.84 


-26.16 


0.21 


0.29 


11 / 144 


800 - 3000 K at OG 


2009/2/6 


2.02 X 0.92 


-16.82 


0.21 


0.17 


15 / 255 




2009/5/18 


2.26 X 0.84 


-23.32 


0.21 


0.28 


15 / 43 




2009/9/13 


2.21 X 0.78 


-19.95 


0.21 


0.25 


- / 54 





" Typical value in self-calibrated images. 

^ The maximum dynamic ranges for phase- referenced (DNp) and self-calibrated image (DN^). 
IR, IS, OG:Iriki, Ishigaki, Ogasawara station, respectively. 



nificantly degraded at this stage. Because of the low 
brightness and large separation angle of J 1755-2232, and 
the low elevation angle during observations, significant at- 
mospheric zenith delay residuals, which is a main cause 
of the coherence loss, still remain. We, therefore, esti- 
mated and corrected this residuals with the image opti- 
mizing method described in Honma et al. (2007, 2008b). 
Estimated residuals were found to be within ±3 cm in 
whole epochs (~ 1.5 cm on average). This is consistent 
with a typical case seen in VERA observations (Honma et 
al. 2008b). The dynamic ranges (or signal to noise ratios) 
of phase-referenced images showed dramatic improvement 
with this correction. But, even after this improvement, 
limited coherence and dcfocusing made several faint spots 
undetectable. 

After these calibrations, synthesized image cube was fi- 
nally made with all maser features in both analyses. Each 
cube had a field of view of 25.6x25.6 mas^ centered on 
relevant feature. Imaging and deconvolution (CLEAN al- 
gorithm) was done in uniform weighting which provided 
the highest spatial resolution. The synthesized beam 
was about 2.2x0.9 mas^ with the position angle of -20°. 
Typical image noise level was ^ 400 mJy beam~^ (I"''') in 
the case of the single-beam analysis (see table 1). It was 
significantly increased in phase-referenced case by a factor 
of 2 or 3, and up to 10 in the worst case. The absolute and 
relative positions of each maser spot were determined by 
an elliptical Gaussian fitting. The formal error in this fit- 
ting was typically 50 /las in RA and 100 /ias in Dec. This 
value is approximately equal to a uncertainty of relative 
positions between each maser spots. Overall discussion 
about the accuracy of absolute positions in each measure- 
ment is given in section 4.2. 

4. Result 

Parallax Measurement 

There were two strong maser features (feature 01 and 
C6 in section 4.3) successively detected for 8 epochs in 
phase-referenced maps. We measured their absolute mo- 
tions referenced on J1755-2232. These motions can be 



expressed by the sum of the annual parallax tt and proper 
motion of each feature /x. The latter is usually assumed 
to be a linear and constant motion for simplicity (e.g., 
Nakagawa et al. 2008). This seems to be applicable for 
our case, because the relative proper motion between two 
features arc actually fitted by a linear motion (see section 
4.3 and figure 2) and the drifts of their line of sight ve- 
locities are well negligible (< 0.2 km s~^ ) during the two 
years. 

A least square fitting was made with tt plus fi to the 
right ascension offsets {X = Aa x cos6 ) from the first 
epoch. We performed our fittings in the same way as 
Hirota et al. (2008). Here, whole detected maser spots in 
each feature (7 maser spots per feature) were distinctly 
used and an initial position of each maser spot, Xq and 
Yq, was included as a fitting parameter. Reduced was 
calculated as 

m — ' 

i 

where m and LUi were the degree of freedom and a fitting 
weight for each data point, respectively. We adopted a 
dynamic range of phase-referenced image (DN) as a fitting 
weight, since it was a rough measure of the coherence loss 
caused by residual delay which was difficult to estimate 
directly. A weight for i-th data point was to be uji oc DN,; 
and scaled to make the total reduced to be ~ 1. If 
we assume the normal distribution, l/y/uJi is equal to the 
error variance of i-th data point. 

Figure 1 shows the examples of the parallax measure- 
ment for G5. 89-0. 39. Note that only a linear proper mo- 
tion was fitted to the declination offsets {Y = AS), since 
G5. 89-0. 39 was located on the ecliptic plane (/3 ~-0.6°), 
where the parallax in the declination degenerated into an 
order of 10 //as. It is just comparable with the maximum 
accuracy of VERA's dual beam astrometry (Honma et al. 
2007) and hard to detect. In each figure, the associated 
error bars indicate l/y^ values. The X were well fit- 
ted by a linear motion (dashed line) and parallax (solid 
curve), by contrast, large dispersion from the best-fit lin- 
ear motion is clearly seen in cases of Y. The error in Y 
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(~^0.9 mas) is, in fact, much larger than that in X 
0.4 mas) as well as several other observations with VERA 
(e.g., Hirota et al. 2007). Detailed analysis about these 
errors is given in next subsection. 

Table 2 presents a summary of these fittings. Presented 
a values are rms deviations of post-fit residuals which 
are good measures of typical error for each measurement. 
Because both of the feature 01 and C6 seem to be asso- 
ciated with a single UCHlI region, we can assume that 
they show the same annual parallax within the astromct- 
ric accuracy of VERA. Thus, the best-fit annual parallax 
was estimated as tt = 0.78 ± 0.04 mas, corresponding a 
distance of 1.28t[! !]^ kpc from the sun, by the weighted- 
means of those for the two features. 

One can see that some data points do not included 
within 1-(T error range in figure 1. This can be due to 
inequable data quality in each observing epoch. We have, 
hence, performed bootstrap analysis to test such an effect, 
where data for each maser spot have been randomly re- 
sampled 10000 times and a least square fitting has been 
iteratively done for each dataset. 

Obtained frequency distributions of tt have been well 
fitted by Gaussian. Detailed parameters are presented in 
table 3. Resultant error of ± 0.05 mas is slightly larger 
than that of simple fittings but almost consistent with 
it, and we finally take it as conclusive l-cr error of our 
parallax measurement. We also emphasize that even if we 
adopt 3-cr error (> 99 % confidence), equivalent distance 
range is still significantly nearer than previously reported 
values. 

4-2. Astrometric Error Sources 

The post-fit residuals seen in our parallax and proper 
motion measurement are much larger than the formal er- 
rors of the elliptical Gaussian fittings. This can be seen 
especially in Y fitting, where residual variance is about 
850 /ias. We again note that the parallax in Y degener- 
ated into negligible order because of the source location. 
This fact suggests that Y fluctuations from a simple lin- 
ear motion directly reflect a degree of actual error in each 
position measurement. One of the main source of this 
error is thought to be atmospheric zenith delay residuals 
as discussed in several studies with VERA (e.g., Sato et 
al. 2008 and references therein). We flrst check this in 
our case using the several expressions in Nakagawa et al. 
(2008). 

They assumed the error in station positions and im- 
age quality of phase-referenced map as additional error 
sources. The former is only about 7 /xas with the separa- 
tion angle of 1.92° and baseline error of 3 mm which is the 
typical value for VERA. This is, of course, negligible in 
our case. The latter can be expressed by a relevant beam 
size in X and Y direction divided by a dynamic range of 
image. The beam position angle are almost parallel to the 
Y direction in our observation, and then, expected error 
is 80 /ias in X and 200 /ias in Y. On the other hand, 
an error caused by atmospheric zenith delay residual is 
proportional to a difference of signal path length between 
2 beams divided by a relevant baseline length. This can 



be calculated from an atmospheric zenith delay residual, 
typical zenith angle and separation angle of the sources 
(see equation (1) in Nakagawa et al. 2008). If we adopt a 
residual of 3 cm as a upper limit in our case and average 
zenith angle of 65°, atmospheric error is about 330 fia.s 
in X and 820 /xas in Y. This is almost consistent with 
the result of Monte Cairo Simulations in Honma et al. 
(2008b). 

Because J1755-2232 looks an ideal point source, the in- 
ternal structure of J1755-2232 does not cause signiflcant 
error. If we also assume that the effect of internal structres 
of maser features is also negligible for simplicity, wc flnally 
gets the total error of 340 and 850 /ias in each direction 
from the root sum square of these 3 values. As one can 
see, this is consistent with the standard deviation of the 
post- fit residuals (see table 2). If differences of a separa- 
tion angle, pair position angle and zenith delay residuals 
are taken into account, this can be roughly comparable to 
the value in Choi et al. (2008), where the source declina- 
tion is almost same in our case. However, the expected 
error of the parallax is estimated to be ~ 0.0 9 mas in a 
statistical way, here, we divide 340 /ias by ^2 x (8- 1) 
considering 8 observing epochs and 2 distinct maser fea- 
tures. This is clearly larger than the actual fitting errors 
of least square analysis (0.04 mas). 

This seems to indicate that the effect of the variability 
of maser feature structures, which we have ignored above, 
cannot be negligible. This type of error is another im- 
portant error source for VERA observations (e.g., Hirota 
et al. 2007). Moreover, it can be actually dominant er- 
ror source in some cases (e.g., Sato et al. 2008). If we 
attribute dominant part of the post-fit variance to this 
effect, instead of the atmospheric error, expected error of 
the parallax can be reduced by the additional factor of 
1/\/N, here is the number of maser spots in single fea- 
ture. Since both of the feature 01 and C6 are consist of 7 
maser spots in our case, the conclusive error of parallax is 
to be 0.33/\/7 x 2 x (8- 1) ~ 0.03 mas, where the numer- 
ator of 0.33 mas is averaged value of the post- fit residual 
in X. This is just comparable to the actual fitting error. 

Therefore, we conclude that the dominant error source 
is the variability of maser feature structures in our case, 
although the contribution from atmospheric zenith delay 
residuals is also nonnegligible. The proportion of their 
contributions is, more quantitatively, to be ^ 2:1, if we use 
averaged atmospheric delay residuals of ~ 1.5 cm which 
has been evaluated from the image optimizing method. 

4.3. Kinematics and Spatial Distribution of Masers 

Absolute proper motions of feature 01 and C6 shown 
in table 2 are motions respect to the Sun, and hence, in- 
clude the contribution of solar motion and the Galactic 
rotation. If we assume the solar motion relative to the 
LSR based on the Hipparcos data (Dehnen & Binney 
1998), the contribution of solar motion is calculated to 
be 0.69 mas yr~^ and -1.19 mas yr""'^ for X and Y, re- 
spectively. The contribution of the Galactic rotation is 
estimated to be -0.04 mas yr^^ and -0.07 mas yr~^ based 
on a i?o of 8.4 kpc and Oq of 254 km s~^ (Reid et al. 
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Feature Ol 



Feature C6 




600 



150 300 450 

Day from 2007/1 1/5 



600 



Fig. 1. The examples of parallax fitting for the brightest maser spot of feature Ol (left) and C6 (right). Here, X and Y axis shows 
elapsed days from the first epoch and positional olfscts from a2000 = 18'^00™30P3066, <52000=— 24°04'04."48649in mas, respectively. 
The annual parallax (solid curve) and linear proper motion (dashed line) were successfully fitted in X direction (upper two panels). 
On the other hand, only a linear proper motion was fitted in Y direction because of the location of G5. 89-0.39 (see text). The errors 
in each panel are calculated from the fitting weights which have been scaled to make the total reduced to be ~ 1. 

Table 2. Summary of simple parallax fittings 



Feature Vlsr (km s ^) 



01 
C6 



TT (mas) 



D (kpc) cr" (mas) (mas yr 



9.4 
11.0 



X 
Y 
X 
Y 



0.79 ± 0.05 
0.77 ± 0.09 



1.27 ± 0.08 



1 30+°-i^ 



0.27 
0.78 
0.39 
0.94 



Combined'^ 0.78 ± 0.04 1.28; 



-0.07 



-0.43 ± 0.09 
-1.43 ± 0.29 
0.77 ± 0.14 
-0.47 ± 0.30 



° The rms deviations of the post-fit residuals. 

^ 1 mas yr~^ — 6.1 km s~^ at the distance. 

" The weighted means of those for two features. 



Table 3. Gaussian Parameters of n distributions from Bootstrap Analysis 



Feature 


Median (mas) 


Standard Deviation 


D (kpc) 


3-f7'' 


01 

C6 


0.80 ± 0.005 
0.75 ± 0.008 


0.06 ± 0.005 
0.10 ± 0.008 


1 25+°-" 

^•^^-0.09 
1 oq+0.20 
^•'^'^-0.15 


+0.38 
-0.24 
+0.87 

-0.38 


Combined'' 


0.78 


0.05 


1 00+0.09 
-'-•'^°-0.08 


+0.33 
-0.22 



" S-fj (> 99 % confidence) error ranges in kpc. 
The weighted means of those for two features. 
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2009) . This value is well negligible and do not change in 
the common range of Rq and 80 (220 - 254 km s~^, 8.0 - 
8.5 kpc;Hou ct al. 2009). Here, wc assumed flat rotation 
and adopted our newly determined distance of 1.28 kpc. 
The intrinsic proper motions of the maser features, where 
these two contributions are subtracted, is represented in 
table 4. We note that these motions are still include the 
peculiar motion of their natal cloud and the internal mo- 
tions of each maser feature. Since only two features are 
detected in 2-beam analysis, we cannot divide these two 
components at this section, but brief presumption will be 
provided in section 5.2.2 based on the model fitting of the 
maser kinematics. 

Additional 12 maser features were detected in the sin- 
gle beam analysis. All parameters of detected features are 
listed in table 5. Total 14 maser features can be divided 
into 4 maser sites. We name these sites as origin (O), 
center (C), north (N) and south (S) based on their posi- 
tions. These four maser sites are widely spread around the 
UCHlI region. We estimated internal proper motions rel- 
ative to the phase- referenced feature 01 for the features 
which were detected in at least three observing epochs. 
Figure 2 shows the example of feature C6. Their internal 
motion was well fitted by linear motion. Estimated mo- 
tions are also summarized in table 4. Overall distributions 
and proper motions of maser features are represented in 
figure 3, where black and grey arrows show relative and 
converted absolute proper motions, respectively. Relative 
positions and internal proper motions of all maser fea- 
tures are converted to absolute values using the position 
and motion of the feature 01. The converted motion of 
C6 is actually coincide with the directly estimated one 
within the error. 

Each of the sites O, N and S, which is located on out- 
side of the ionized shell, contains only one or two maser 
features. This limited number of maser features can be 
related to a life time of H2O maser activity. They are 
generally believed to disappear along with a evolution of 
UCHlI region (e.g., Beuther et al. 2002b, Breen et al. 

2010) , in attributing to a dispersion of dense gas which 
required to excite masers. All of these sites are associated 
with SiO J = 8 - 7 emission detected by Sub-Millimeter 
Array (SMA) (Hunter et al. 2008). SiO J = 8 - 7 emis- 
sion is a moderate shock tracer and frequently observed 
in post-shock gas which is associated with a protostellar 
outflow (e.g., Takami et al. 2006). Such a spatial relation- 
ship may indicate that H2O maser and SiO emission trace 
same shock fronts. This is not so surprising case because 
H2O masers are also associated with strong outflow shock 
frequently (e.g., Motogi et al. 2008; Torrefies et al. 2010). 

On the other hands, line of sight velocities of maser 
features in the site O and N are not consistent with that 
of SiO emission. The velocity offsets are larger than 5 
km s^^. These offsets may reflect the difference of pre- 
cise locations where each emission comes from. This can 
be simply attributed to the different excitation conditions 
of SiO and H2O. The critical density of SiO J 8 - 7 
emission are tihj — 10^ cm~^ based on the database in 
Schoier et al. (2005). It is actually 2 orders of magnitude 



smaller than that of H2O maser and comparable with pre- 
shock density of the maser (Elitzur 1992). H2O maser is, 
hence, probably excited in the most strongly compressed 
part such as a head part of a bow shock. Emission from 
that region cannot dominate the integrated SiO emission, 
since such a region should have quite limited volume com- 
pared to total post-shock gas. 

The site C, by contrast to other site, is not associated 
with SiO emission. There are 9 maser features and almost 
all of them (7 of 9) are highly variable and only detected in 
a single epoch. This maser site had been reported in the 
past VLA observation (Hofner & ChurchwcU 1996), but 
observed position and line of sight velocity was slightly 
different from our detection. Whole maser features, in- 
cluding the one detected by VLA, closely located at the 
position of Feldt's star (< 200 mas). Their line of sight 
velocities are little red-shifted from systemic velocity of 
9 km s~^ (e.g.. Hunter et al. 2008) and clearly different 
from that of OH masers seen in the same position (Stark 
et al. 2007). 

Highly blue shifted (~ -35 km s^^) OH masers are 
thought to be excited in expanding neutral shell or strong 
outflow (the former in Stark et al. 2007, the latter in 
Zijlstra et al. 1990). They seem to be located on the 
foreground of the ionized shell in either cases. The kine- 
matic difference strongly invokes the distinct origins of 
OH and H2O masers. Such a nearly located, but strictly 
distinct displacement of these two maser was quite natural 
if we considered different excitation conditions again, and 
it had been actually observed for several regions (Forster 
& Caswell 1989). Taking into account these situation, we 
propose that this H2O maser site is just located inside a 
ionized shell and excited in a remnant of dense circumstel- 
lar structure such as an accretion disk. Extremely young 
dynamical age of UCHlI region 600 yr) well support 
the presence of such a remnant and their crescent-like dis- 
tribution and velocity field can be actually explained by 
a simple ring model (see section 5.2.2). 

5. Discussion 

5.1. Validity of New Distance 

The newly estimated distance of 1.28 kpc is well smaller 
than the previously reported distances (> 1.9 kpc). This 
is mainly because almost all of past measurements have 
been done kinematically. Kinematic distance essentially 
depends on the Galactic rotation model and accuracy of 
determining rotation velocity of the source. If it is the case 
that target source has unknown, non-negligible velocity 
component aside from the Galactic rotation, kinematic 
distance is no longer reliable in principle. This type of 
systematic error is more and more significant especially 
for the source near the Galactic center direction. In the 
case of our source {I = 5°. 89), only a few km s^^ offset 
along the line of sight can cause an error of ~ 1 kpc. 

Rcid ct al. (2009) had first analyzed the Galactic ro- 
tation based on the distances of 18 star-forming regions 
which are estimated only from trigonometric parallaxes 
of masers derived from VLBA and VERA observations. 
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Table 4. Intrinsic Proper Motions 

Feature Absolute Proper Motion" Relative Proper Motion' 



fJ-x tJ-Y fJ-x fJ-Y 



01" 


-1.08 ± 0.24 


-0.16 ± 0.40 






C6 


0.12 ± 0.29 


0.79 ± 0.47 






02" 


-1.02 ± 0.33 


-0.46 ± 0.43 


0.06 ± 0.09 


-0.29 ± 0.03 


C6 


-0.21 ± 0.29 


1.31 ± 0.48 


0.87 ± 0.05 


1.47 ± 0.08 


C9 


-1.51 ± 0.27 


0.72 ± 0.53 


-0.43 ± 0.03 


0.89 ± 0.13 


N 


-3.38 ± 0.43 


-1.79 ± 1.56 


-2.3 ± 0.19 


-1.63 ± 1.16 


SI 


-0.12 ± 0.40 


-1.03 ± 0.93 


0.96 ± 0.16 


-0.86 ± 0.53 


° All values are in the units of mas yr 


-i 





The upper 2 rows are from the parallax fitting. 

The lower 5 rows are from the internal proper motions relative to 01. 




Fig. 3. Overall distribution of mascr features. Lower left panel: The locations of 4 maser site (0,C,N,S) are presented by color points 
with dashed circle. The coordinate origin is a2000=18'^00™30?3, i52ooo=~24°04'0."0. Here, background color scale, black crosses 
and contour shows the first moment map of SiO J = 8-7 emission, 875 fj,m dust cores and 3 mm continuum image (0.1 Jy beam~^) 
from Hunter et al. (2008), respectively. The star marks the positions of Feldt's star a2000=18'"00™30?44, 52OO0=-24°04'0f'9. Other 

4 panels: Detailed distributions of maser features are shown. Here, each point indicates detected maser feature. Each axis shows 
the relative coordinate from the phase-referenced maser feature Ol in the units of mas. Absolute and internal proper motions are 
written in black and grey arrows, respectively. Red triangle with error bar in site C (upper right panel) is the maser detected by 
Hofncr & Churchwell (1996). Here, the star marks the position of Feldt's star again and its error is presented by dashed circle of 
250 mas diameter (200 from positional error and 50 from PSF, see Feldt et al. 2003). The color represents Vlsr (km s~^) for all of 

5 panels in the scale of associated color bars. 
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Table 5. Detected Maser Features 



Feature 


Epoch 


^SR 




Offset" (mas) 




Peak Intensity 
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5 
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3.9 


N 


2 


3.65 


1.05 
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6.7 




3 


2.95 


2.95 


3932.14 (0.02) 


4186.54 (0.04) 


18.0 




7 


3.82 


1.26 


3929.04 (0.05) 


4184.55 (0.11) 


2.9 


SI 


4 


8.23 


1.26 


4898.62 (0.07) 


-1070.77 (0.19) 


8.0 




5 


8.11 


1.05 


4898.48 (0.04) 


-1070.66 (0.10) 


13.7 




6 


8.06 


1.47 


4898.27 (0.05) 


-1071.59 (0.15) 


12.fi 




7 


8.03 


1.47 


4898.35 (0.03) 


-1071.68 (0.07) 


11.3 




8 


8.13 


1.05 


4898.50 (0.05) 


-1071.14 (0.14) 


4.7 


S2 


8 


8.13 


0.63 


4897.28 (0.09) 


-1075.02 (0.24) 


2.4 



° Full width at zero intensity (FWZI) for each maser feature. 

" The positions relative to (18''00'"30f30fifi, -24°04'4"48fi49) (J2000.0). 
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Fig. 2. The internal proper motion of feature C6. Y axes in 
each panel show the positional offset relative to the phase-ref- 
erenced feature Ol. Both of Dashed lines indicate the best-fit 
linear proper motions. 

Their new analysis indicates that commonly used kine- 
matic distances are generally overestimated, sometimes by 
factors greater than 2, which is just in our case. Similar 
result has already been reported in other VERA obser- 
vation (G14.33-0.64, the measured distance ^ 1.1 kpc, / 
= 14°. 33 (Sato et al. 2010a)). Our result seems to be 
well consistent with the hypothesis in Sato et al. (2010a) 
that the Sagittarius spiral arm, in which G5. 89-0. 39 and 
G14.33-0.64 are thought to be located, lies at the closer 
distance (^ 1 kpc) from the sun compared to the previous 
value (2-3 kpc). Of course, more large samples must be 
required to confirm it. 

Progressive Galactic simulation in Baba et al. (2009) 
also suggests that star-forming regions and young stars in 
a spiral arm have significant non-circular motions up to 
30 km s^^. If their calculation is correct, a kinematic dis- 
tance intrinsically has large systematic error of 2-3 kpc. 
With these contexts, the significant difference between the 
kinematic distance and our new distance can be attributed 
to the systematic error in the kinematic distance. 

ACW98 have been derived the source distance of 2.0 kpc 
from the expansion of the ionized shell, but it is still highly 
model dependent (see section 5.2.3). Consequently, our 
new distance which has no model dependency is thought 



to be the most reliable distance of G5. 89-0. 39. 

5.2. Physical Parameters and High Mass Star Formation 
m G5. 89-0.39 

Because of significant modification of the distance 
to G5. 89-0. 39, we first recalculated physical parameters 
which had been previously reported based on the new dis- 
tance. Table 6 to 9 summarize recalculated parameters. 
Each table contains both of original and modified values 
with the reference papers. Detailed explanations for listed 
parameters are given in related subsections and captions. 
5.2.1. Evolution of the 0-type Protostar 

Table 6 show the properties of Feldt's star. We de- 
rived stellar spectral type in the same manner in Wood & 
Churchwell (1989). We first estimated excitation param- 
eter of the source from radio continuum emission (see also 
Krutz et al. 1994). This was simply converted to the flux 
of Lyman continuum photons (A^l) and spectral type (Sp 
type) comparing the stellar models in Panagia (1973). In 
the second, we used total far-infrared luminosity (ipm)- 
LpiR can be a good measure of bolometric luminosity for 
deeply embedded sources. Because of the low spatial reso- 
lution, LpiR may include several contributions from other 
objects around the UCHlI region, and hence, this gives 
upper limit of the source luminosity. 

Derived spectral types arc almost consistent and indi- 
cate late 0-type star (08 - 08.5 ZAMS). Expected steUar 
mass A/* is about ~25 Mq (e.g., Vacca et al. 1996). This 
is a roughly half value of previously reported 05 ZAMS 
(or 05V) star. Feldt et al. (2003) tested such an early 
spectral type by a model fitting of mid-infrared color {Ks 
and L' bands). But their fitting is based on the distance 
of 1.9 kpc, and if we adopt new distance and appropriate 
extinction value, all of O type stellar models which are 
later than 05 can move into the range of L' band excess 
reported in Feldt et al. (2003). 

The line ratio of He I triplet at 2.11 fim to Br7 also sup- 
ports the existence of late 0-type star. Hanson, Luhman, 
& Rieke (2002) have observationally determined upper 
limit of this ratio as 0.02 for G5. 89-0. 39 and indicated 
that this value is clearly smaller than a value expected 
for a star earlier than 07 type. Puga et al. (2006) actu- 
ally detected 2.11 ^m He I emission with VLT. Although 
they suggested rather early O star (< 07), the detection 
only at the location of Feldt's star also ruled out a star 
earlier than 07, because such an early type star could 
fully ionize whole He atoms and provide constant line ra- 
tio throughout a UCHlI region (e.g., Ostcrbrock 1989). 
This seems to be a reliable constraint, since the line ratio 
is fully independent of the source distance. Gonsequently, 
new spectral type of 08 - 08.5 seems to be reasonable. 

As seen in table 6, we also estimated the momentum 
rate of stellar wind (Avind)- It has been calculated from 
pressure balance at the inner surface of ionized shell. We 
simply assumed that spherically-symmetric wind and then 
thermal pressure of ionized gas was balanced with Pwind 
per unit area. Although detailed properties of mass loss 
activity at the ZAMS stage are still unknown, estimated 
-Pwind is ten times larger than that of a dwarf star which 
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Table 6. Recalculated Stellar Parameters 



Parameter 


Original value 


Modified value 


Reference 




4.5xl0''« 


IJxlO-** 


(1) 


^FIR {Lq) 


3.0x10^ 


7.3x10-* 


(1) 


Sp Type (radio) 


07 


08.5 


(1), (2) 


Sp Type (FIR) 


06 


08 


(1), (2) 


(Mq) 


~ 40 


-25 


(3) 


Pwind (Mekm s-i yr-1 ) 




2.5 X IQ-^ 


(4) 



References.- (l):Wood & Churchwell 1989, (2):Panagia 1973, 
(3):Vacca et al. 1996, (4):this work. 



Table 7. Recalculated Parameters of UCHll region 



Parameter 


Original value 


Modified value 


Reference" 


Tie (cm^"^) 


2.4x10^ 


3.2x10^ 


(1) 


i?i/-Rc (au) 


1900 / 5200 


1200 / 3300 


(2) 


Ushcu'' (km s"^) 


39 


25 


(2) 


tdyn (yr) 


600 


600 


(2) 


-Bxherm (10^^ Org) 




0.7 


(3) 


EKin (10^5 erg) 




1.2 


(3) 


Pshcii (Mgkm ) 




5.6 


(3) 


Monv (Mo) 


300 


123 


(4) 


rienv (cm"3) 


5.3 xlO^ 


8.0 xlO^ 


(4) 


Tenv (K) 


40 - 140 


40 - 140 


(5) 



"References.- (l);Wood & Churchwell 1989, (2):ACW98, (3):this work, 
(4):Tang et al. 2009, (5):Su et al. 2009. 
Tangential expansion derived from the angular expansion rate. 



Table 8. Parameters of the E-W Outflow (CO J =1 - 0)) from Watson et al. (2007) 



Parameter 


Original value 


Modified value 


tdyn (yr) 


7700 


5000 




Blue lobe Red lobe 


Blue lobe 


Red lobe 


"flow (km s-^) 


-15.0 - 3.2 13.2 - 25.0 


-15.0 - 3.2 


13.2 - 25.0 


M (Mo) 


123 116 


50.4 


47.5 


-Bkin (10^6 erg) 

P (10^ Mokm s-i) 


6.4 5.3 


2.6 


2.2 


9.0 8.0 


3.7 


3.3 


M (10-3 Moyr^i) 


16.0 15.0 


10.0 


9.6 


P (10~2 Afgkm s-i yr-i ) 


12.0 10.0 


7.7 


6.4 


-^mcch (-^0) 


68.0 56.0 


44.0 


36.0 



" Integrated velocity range. 



Table 9. Parameters of the inner high velocity Outflow (CO J =3 - 2)) from Klaassen et al. (2006) 



Parameter 


Orig 


inal value" 


Modified value 


^dyn (yr) 


2000 


1300 




Blue lobe Red lobe"^ 


Blue lobe 


Red lobe'^ 


"flow (km s-^) 


-66.0 


78.0 


-66.0 


78.0 


M (Mo) 


2.7 


0.6 


1.1 


0.25 


i;ki„ (10^6 erg) 


1.5 


0.3 


0.6 


0.1 


P (Mokm s-i) 


55.9 


12.0 


22.0 


4.9 


M (10-4 Moyr-i) 


13.5 


3.0 


8.5 


1.9 


P (10^2 Mokm s-i yr-i ) 


2.8 


0.6 


2.7 


0.4 


imcch (-^0) 


59.0 


11.5 


38.0 


7.4 



" Inclination of outflow axis which has been originally assumed as 45° is ignored 
The maximum velocities. 

Absorption features seen in the red lobe have not been corrected 
in Klaassen et al. (2006). 
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has ^25 A/0 (Sternberg et al. 2003 and reference therein). 
A wind velocity is generally determined by gravitational 
force at the launching point of a wind, and it is thought 
to be same order between ZAMS and dwarf stage. This 
is about 2.0 xlO^ km s"^ for 25 Mg dwarf (e.g., Smith 
et al. 2002). In this case, mass loss rate is order of 10~^ 
M© yr~^, and of cource, ten times larger than the value 
of dwarfs . One possible explanation for this excess may 
be larger stellar radius of very young ZAMS star which is 
still under contraction. This could cause relatively weak 
surface gravity and might help mass loss activity. 

Recent theoretical works in Hosokawa & Omukai 
(2009) calculated detailed stellar evolution via spherically- 
symmetric mass accretion. They have shown that large 
entropy supply under high accretion rates (> 10~^ Mq 
yr~^) onto a proto-stellar surface prevent contraction of 
protostar and, at some stage, cause far large stellar radius 
(- 100 Rq). Hosokawa, Yorke, & Omukai (2010) has 
been treated the case of cold-disk accretion, and they have 
suggested that expected evolution is qualitatively same as 
spherically-symmetric case after 10 M©. Realistic mass 
accretion should correspond to some intermediate state 
between these two extreme case as they pointed out. 

The most important product of their evolutional track 
is delayed hydrogen burning. Resultant young high mass 
(proto-)star with large radius and limited UV photons 
can explain highly luminous but low-effective temperature 
source which has no detectable Hll region (e.g., Beuther 
et al. 2007, Furuya & Shinnaga 2009). If Fcldt's star 
has been just after the hydrogen ignition, the maximum 
accretion rate, with which the star has been evolved, is 
estimated to be ^ 10^'^ Af© yr^^ from ^ 25M© based 
on their calculations. This assumption may mean that 
Fcldt's star is quite rare sample, but it seems to be well 
consistent with extraordinarily young dynamical age of 
UCHlI region. If this is the case and we assume steady 
mass accretion, the time required to form Fcldt's star is 
simply ~ 2.5 x 10"* yr. 

5.2.2. Maser Morphology Near the Feldt's Star: a 
Possible Partial Ring 

Newly detected crescent-like structure of water maser 
may also suggest youth of the source. As we mentioned 
above, the maser clumps in the site C can be associated 
with circumstellar remnant gas. Their alignment and ve- 
locity field are actually fitted by the Keplerian rotating 
and expanding ring model in Uscanga et al. (2008). The 
fitting parameters are central position of a ring (xo,yo), 
radius i?, inclination from the celestial plane i, position 
angle of apparent major axis measured from east toward 
north 6, expansion (or infalling) velocity Vexp and velocity 
of the source {v^, Vy, Vz). Here, z means a line of sight 
direction. We regarded Feldt's star as a central source 
of ring, in turn, its mass was fixed to be 25Af0. Thus, 
Keplerian rotation velocity depended on R only. 

We performed the model fit with a standard fitting. 
First, the spatial distribution of 9 (our detection) -I- 1 
(VLA detection) maser features was fitted by xq, yo- R- h 
6, and then, v^, Vy, were determined from the velocity 
field of masers. It should be noted that and Vy are 



Table 10. The best fit parameters of the ring model 



Parameter 


Best Fit (err°) 


xQ (mas) 


1914.6 (0.5) 


yO (mas) 


3877.2 (0.8) 


R (mas) 


235.2 (1.2) 


en 


110.5 (0.2) 


i" (°) 


+83.0 (0.2) 


•Uexp (km s ^) 


-1.2 (0.8) 


(km s~-^) 


-6.5 (1.5) 


Vy (km s""'^) 


6.4 (1.1) 


(km s^^) 


18.6 (0.4) 



" Formal 4 a error of fittings. 
Western edge is on the far side. 



more uncertain than , because there are only 2 available 
proper motion vectors. Table 10 and figure 4 show the 
best fitted parameters and estimated ring superposed on 
the maser alignment, respectively. 

The best fitted ring is almost edge-on and its position 
angle is roughly north - sorth direction. It is noteworthy 
that this geometry matches well the largest outfiow which 
has a position angle of nearly linear alignment from east 
(red lobe) to west (blue lobe) (see below). The center 
of the ring is slightly offset from Feldt's star (1826 mas, 
3587 mas), but it can be within the ring if we take into ac- 
count ~ 200 mas error. The small and negative expansion 
velocity means that the ring is slowly infalling. 

Class II CH3OH masers which are believed to trace an 
accretion disk around HMPOs (e.g., Bartkiewicz et al. 
2009 and reference therein), are not seen in G5. 89-0. 39 
region. It may be consistent with the fact that Class 
II CH3OH maser emission trace pre-UCHlI stage (e.g., 
Walsh et al. 1998; Codella & Moscadelli 2000). On the 
other hand, a disk-like structure traced by 22 GHz H2O 
masers has been actually reported in some sources. Imai 
et al. (2006) has especially found infalling-rotating maser 
disk, which is similar case to ours, for early B-type HMPO. 
Torrelles et al. (1996) have reported the H2O maser disk 
candidate in Cepheus A HW2, which had almost same ra- 
dius (^ 300 au) as our case. If the crescent-like masers 
actually traces an accretion disk (or its remnant), it is 
clear indicative of a quite youthful 0-typc object which is 
at least before complete disk evaporation. 

We finally discuss about the velocity field. The 3D mo- 
tion of the ring [vx, Vy, Vz) still includes the non-circular 
(peculiar) motion of G5. 89-0. 39 system (recall section 4.3 
and 5.1). To divide these contributions, we now assume 
2 conditions as follows, (1) Feldt's star has been initially 
located at the center of the radio shell and migrated to 
the current position in the dynamical time of the shell, 
(2) the ring moved together with the star. The possibility 
of stellar migration was already discussed in Feldt et al. 
(2003). This motion corresponds to be a proper motion of 
about 1.04 ± 0.30 and 0.98 ± 0.28 mas yr"! in X and Y 
direction, respectively. By contrast, the migration veloc- 
ity along a line of sight (~ 9.6 km s~^) can be obtained 
from Vz - VhSK- The 3D motion of the partial ring can 
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Fig. 4. Left panel shows the best-fit ring model superposed on the maser distribution in figure 3. Curved arrow represents rotational 
direction. The red and blue arrows show a same geometry as red and blue lobe of E-W outflow. The inclination of the ring is consistent 
with this outflow geometry. Right panel shows velocity field of the ring. Red and black arrows show expected proper motion vectors 
of the ring and observed vectors, respectively. These motions contain both of the systemic velocity of G5. 89-0. 39 offset from the 
Galactic rotation and migration velocity of the ring. 



Table 11. Divided velocity components 



Migration velocities of Feldt's star 





6.3° (0.3) 6.0 (0.3) 9.6 (0.0) 


Non-cireular motion of G5. 89-0. 39 


(f/, y, W) 


-12.8 (1.8) 0.4 (1.4) 4.7 (0.4) 
5.4 (1.6) -5.5 (2.4) 11.3 (2.4) 



Velocities are in the units of km s 



be eventually divided into the migration motion of Feldt's 
star and systemic peculiar motion of G5. 89-0. 39. These 
divided values are presented in table 11, where the line of 
sight component of the peculiar motion (~ 4.7 km s~^) 
was also extracted from Vlsr using the Galactic rotation 
model in Hou et al. (2009). 

The peculiar motion of the system have also been shown 
in ([/, V , W) frame at the same time. Here, U is toward 
the Galactic center from the source position, V is along the 
Galactic rotation at the local place, W is perpendicular 
to the Galactic plane and toward the north Galactic pole. 
Magnitude of the non-circular component ~ 14 km s~^ is 
enough within theoretically predicted range (see section 
5.1 again). 



5.2.3. Anisotropic Expansion of UC'HlI Region 

Table 7 summarizes the parameters of the UCHlI re- 
gion. Upper seven rows present the electron density rig, 
inner / outer radius of the radio shell i?i / Ro, expansion 
velocity Wshcii, dynamical age tdyn, thermal energy of ion- 
ized gas -Ethcrm, kiuctic energy i?kin and total momentum 
of shell expansion Pshcii. Lower three rows contain the pa- 
rameters of surrounding dense envelope. Envelope mass 
Mcnv and density obtained from SMA observation 

of submillimeter dust continuum (Tang et al. 2009) and 
its size scale is roughly 2 x 10^ au. In especial, envelope 
temperature Tonv have been derived from an excitation 
analysis of the CH3CN (J = 12 - 11) line by Su et al. 
(2009), and hence, do not depend on a distance. They 
have found a temperature gradient in the range shown in 
the table. 

Overall properties of the ionized region are not changed 
drastically except for the expansion velocity which is cal- 
culated from angular expansion rate of 4 ± 1 mas yr~^ in 
ACW98. Even if ionized gas is in a free expansion state, 
its expansion velocity should be almost equal to a sound 
speed of ionized gas (^ 13 km s^^ at 10* K). This sug- 
gests that only a half of total momentum can be thermally 
provided and the shell must require another source of mo- 
mentum input. If we attribute it to a spherically symmet- 
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ric stellar wind discussed in section 5.2.1, required time 
is about 1000 yr. Moreover if the momentum of photons 
from Feldt's star is added, the total momentum of the 
shell could be provided during idyn- 

ACW98 have determined the source distance in compar- 
ison between angular expansion rate and the line width 
of radio recombination lines under several assumptions 
about how the shell expansion contributes to the line 
width. Now we have more accurate distance, and in turn, 
we can estimate a residual contribution to the line width 
which offsets from spherical expansion. This value is 
50 km based on their analysis. This is twice as large 
as tangential expansion and nearly consistent with highly 
blue shifted OH masers in front of the shell. 

This highly anisotropic expansion could be explained 
by the idea proposed by Hunter et al. (2008), where an 
outflow from SMAl, which is the luminous dust core near 
Feldt's star, had been going across north to south and it 
was disrupted after the formation of the UGH II region. If 
this is the case, ionization of intrinsically expanding gas 
within the outflow can give a plausible explanation for 
high velocity ionized gas, and hence, observed anisotropy. 
5.2.4. Outflow Nature 

The revised properties of the outflowing gas in G5.89- 
0.39 are listed in tables 8 and 9. Table 8 shows the pa- 
rameters ofGO(J = l- 0) emission which trace the outer 
extended flow (Watson et al. 2007), and table 9 presents 
that of the inner high velocity flow which has been de- 
tected in CO ( J = 3 - 2) emission (Klaassen et al. 2006). 
Each table contains the dynamical age tdym outflow ve- 
locity (or velocity ranges) vrow, outflow mass M, kinetic 
energy £^kin, momentum P, outflow rate M, momentum 
rate P and mechanical luminosity imoch- We note that 
an effect of inclination is not corrected in these tables, 
although Klaassen et al. (2006) have originally assumed 
inclination of 45°. 

The largest outflow in table 8 is thought to be driven 
by Feldt's star and extends just east - west direction (see 
figure 2 in Watson et al. 2007). Total mass measured from 
CO (J = 1 - 0) is ~ 100 A/q and greatly exceeds the 
source mass. Such an extremely massive outflow is often 
observed toward a high mass star-forming region (e.g., 
Lopez-Sepulcre et al. 2009). The origin of large mass is 
unsolved problem (e.g., Churchwell 1997). 

A detailed driving mechanism of a protostellar outflow 
is still under discussion, however Machida, Inutsuka, & 
Matsumoto (2008) have shown that high velocity jet and 
molecular outflow (hereafter called intrinsic outflows) are 
distinctly driven in their MHD simulation as natural prod- 
ucts of magnetized core collapse and mass accretion via 
rotating disk. This is consistent very well with a case of 
nearby low mass star formation, and seems to be appli- 
cable for high mass case as long as objects are formed 
by gravitational collapse and mass accretion. The energy 
source of intrinsic outflows is gravitational energy which 
is transformed into outflow energy mediated by magnetic 
fields. In this point of view, total outflow mass of intrinsic 
outflows cannot exceed that of central object at a driving 
point. Clearly, much of the outfiow mass arises from mate- 



rial entrained far from the driving source. Quantitatively, 
observations of outfiows from massive young protostellar 
objects show that the entrained mass is on the order of 
4% of the core mass surrounding the central protostar 
(e.g., Beuther et al. 2002a). This relation suggests that 
the parent core for the forming cluster containing Feldt's 
star contained about ~ 2500 Mq. 

On the other hand, if we assume the ratio of the outflow 
to the accretion rate as '--^ 10 % from theoretical works (e.g, 
PcUetier & Pudritz 1992), total mass of intrinsic outflows 
for Feldt's star is ~ 2.5 Mq following the discussion in 
section 5.2.1. This is comparable with the mass of the 
GO(J = 3- 2) outflow which traces inner hot and high 
velocity outflow (see table 9). If we consider a momentum 
transport from the GO (J = 3 - 2) outflow to the GO (J 
= 1-0) outflow, required time is 2.3 x 10^ yr. This is 
consistent with a formation time of Feldt's star estimated 
above and clearly larger than the dynamical age of outflow 
which is obtained from an extent divided by a current 
velocity. We note that this time scale should be a lower 
limit, since the mass and momentum of the CO ( J = 3 - 
2) outflow is the sum of the east - west flow from Feldt's 
star and north - south flow from SMAl. They have been 
unresloved in the single dish data in Klaassen et al. (2006), 
and actually, significant fraction of high velocity emission 
arises from the latter (see high resolution images in Hunter 
et al. 2008). 

The fact that a commonly used dynamical time under- 
estimates actual outflow age has already been pointed out 
in the case of low mass objects (Parker et al. 1991) and we 
suggest that it can occur in the case of too massive outflow 
in high mass star formation which includes large entrained 
mass. It should be, consequently, careful that an observa- 
tionally estimated outflow rate, which often reaches 10~^ 
Mq yr~^ for high mass star, is also overestimated by 
an order of magnitude. This looks reasonable because 
Hosokawa & Omukai (2009) have predicted that large ac- 
cretion luminosity must prevent steady accretion in the 
case of such a too large accretion rate. Detailed infor- 
mation about a termination of high velocity inner flow is 
quite useful to conflrm these discussions. This would be 
obtained from direct measurement of proper motions of 
outflow lobes with ALMA. 

6. Conclusion 

The distance of G5. 89-0. 39 is newly estimated to be 
1.28lQ Qg kpc from the annual parallax measurement with 
VERA. This is 2/3 of the previously known value, but it 
is well reasonable if we take into account the small galac- 
tic longitude of ^ 5.89° and recent theoretical prediction 
about non-circular motions of star-forming regions. 

Rescaled physical parameters based on the new distance 
give us several in-depth natures of high mass star forma- 
tion in G5. 89-0. 39 as follows. 

(1) The ionizing star is rather later type ZAMS than 
previously believed type of 05. Spectral type of 08.5 - 
08 means that the UCHlI region arc excited by not so 
massive and standard 0-type object. Expected accretion 
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rate is ~ 10^^ Mq yr~^ based on the extremely young 
age of the ionized shell and detailed evolutional track of 
massive protostar under a high accretion rate. Resultant 
formation time is about ^ 2.5 x IC^ yr in this case. 

(2) Detected maser alignment at the 0-star can be fit- 
ted by infalling Keplarian ring and its inclination and po- 
sition angle are also consistent with east-west orientation 
of the strong outflow. It seems to be trace accretion disk 
(or its remnant) and suggest remarkable youth of the O- 
star which is before complete evaporation of circumstellar 
structure. 

(3) Reconsideration of outflow nature suggests that 
the large portion of outflow mass should be entrained 
from massive envelope. A commonly used dynamical 
time should significantly underestimate actual outfiow age 
same as low mass cases (Parker et al. 1991). This also 
causes an overestimate of outfiow rate by an order of mag- 
nitude. Direct observation of momentum transportation 
from intrinsic outfiow to outer entrained flow is required 
to confirm this. This may be able to be achieved with the 
proper motion measurement of outflow lobe with ALMA. 

We finally emphasize that G5. 89-0. 39 is one of the 
nearest target to investigate individual high mass star 
formation and evolution of core scale cluster including an 
0-type object. 
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